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Abstract. Thermal models are commonly used to interpret heavy- ion data on 
particle yields and spectra and to extract the conditions of chemical and thermal freeze- 
04 | out in heavy- ion collisions. I discuss the usefulness and limitations of such thermal 

model analyses and review the experimental and theoretical evidence for thermalization 
. in nuclear collisions. The crucial role of correlating strangeness production data with 

^ ' single particle spectra and two-particle correlation measurements is pointed out. A 

consistent dynamical picture for the heavy-ion data from the CERN SPS involves 
| an initial prehadronic stage with deconfined color and with an appreciable isotropic 

■ pressure component. This requires an early onset of thermalization. 
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O ' 1. Introduction 
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The stated goal of the relativistic heavy-ion programs at CERN and BNL is the study of 
the phase diagram of strongly interacting matter at high temperatures and densities and 
the search for the quark-gluon plasma (QGP). The discussion of a phase diagram requires 
thermodynamic language. A phase transition from an initial color-deconfmed QGP to a 
color-confined hadronic state (as it is supposed to occur in heavy-ion collisions) can only 
be reasonably well defined if the system under study is in a state of approximate local 
thermodynamic equilibrium. The application of thermal and hydrodynamic models to 
relativistic heavy-ion data is therefore more than a poor man's approach to heavy-ion 
dynamics, it is rather a necessity for everybody who wants to convince himself and 
others that we succeeded in creating the quark-gluon plasma (QGP) and observing the 
phase transition accompanying its hadronization. 

Of course, the models may fail; in fact, they must necessarily fail beyond a certain 
level of detail when applied to heavy-ion data. The reason is obvious: the collision 
systems are small, causing corrections to the infinite volume limit usually assumed in 
the thermodynamic approach, and they undergo a strong dynamical evolution on time 
scales which are comparable to the microscopic thermalization time. Thermal models 
therefore can never provide more than a rough picture of the bulk of the phenomena, 
good for qualitative answers; on a more detailed and quantitative level, the failure 
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of the thermal model will become manifest, and traces of genuine QCD dynamics (as 
opposed to thermodynamics) will show up. But when trying to assess bulk phenomena 
like QGP formation, we are not (in fact, we must not be) primarily interested in these 
deviations from thermodynamic behaviour and the traces of elementary QCD dynamics; 
the latter can be studied much more easily and cleanly in elementary lepton or hadron 
collisions. We should rather concentrate on the rough global pattern of the data and try 
to understand them within a (nota bene: sufficiently sophisticated, see below) thermo- 
and hydrodynamic approach. On the other hand, if it turns out that not even the rough 
qualitative features of the data can be understood in this way and detailed hadronic 
dynamics is required even for a superficial understanding of the observations, then we 
should concede that our attempt to create "hot QCD matter" has failed. 

Having argued in favor of a "simple" thermal approach to heavy-ion data, the next 
questions to be addressed are (i) the level of sophistication which the thermal model 
should have before being applicable to the description of particle production in nuclear 
collisions, (ii) which level of agreement between model and data can at best be expected, 
and (iii) where to draw the line between agreement and disagreement when comparing 
model and data. This is what this contribution is about. At the end I will try to 
draw some conclusions about what we have learnt so far from the thermal analysis of 
heavy-ion data, and which further steps should be taken. 

2. Two types of "thermal" behaviour 

"Thermal" behaviour can arise in conceptually different ways, with different meanings of 
the "temperature" parameter T. For us the two most important variants of "thermal" 
behaviour are the following: 

1. The statistical occupation of hadronic phase-space with minimum information. 
The latter is in practice provided by external constraints on the total available energy 
E, baryon number B, strangeness S and, possibly, a constraint A s on the overall fraction 
of strange hadrons. "Thermal" behaviour arises in this case via the Maximum Entropy 
Principle in which the "temperature" T and "fugacities" e w " /T , e^ 3 ^ (which in the 
canonical approach are replaced by so-called "chemical factors" |], occur as Lagrange 
multipliers for the constraints. Examples are nucleon emission from an evaporating 
compound nucleus in low-energy nuclear physics and hadronization in e + e~, pp and pp 
collisions (hadron yields ffl, [| and m_|_-spectra ||). The number of parameters to fit the 
data in such a situation is equal to the number of "conserved quantities" (constraints), 
and it reflects directly the information content of the fitted observable (s). This type of 
"thermal" behaviour requires no rescattering and no interactions among the hadrons, 
there is no isotropic pressure and no collective flow in the hadronic final state and, in fact, 
the concept of local equilibrium can not be applied. Of course, this type of "thermal" 
behaviour is not really what we are interested in in heavy ion collisions, except as a 
baseline against which to differentiate interesting phenomena. 
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2. Thermalization of a non-equilibrium initial state by kinetic equilibration 
(rescattering) . This does require (strong!) interactions among the hadrons. Here 
one must differentiate between thermal equilibration (reflected in the shape of the 
momentum spectra), which defines the temperature T, and chemical equilibration 
(reflected in the particle yields and ratios) which defines the chemical potentials in 
a grand canonical description. The first is driven by the total hadron-hadron cross 
section while the second relies on usually much smaller inelastic cross sections and thus 
happens more slowly. This type of equilibrium is accompanied by pressure which drives 
collective flow (radial expansion into the vacuum as well as directed and elliptic flow in 
non-central collisions). In heavy ion collisions it is realized at most locally, in the form 
of local thermal and/or chemical equilibrium - due to the absence of confining walls 
there is never a global equilibrium. This is the type of "thermal" behaviour which we 
are searching for in heavy-ion collisions. 

I stress that flow is an unavoidable consequence of this type of equilibration. 
Thermal fits without flow to hadron spectra are not consistent with the kinetic 
thermalization hypothesis. Flow contains information; it is described by three additional 
fit parameters v(x). This information is related to the pressure history in the early stages 
of the collision and thereby (somewhat indirectly) to the equation of state of the hot 
matter. 

Most thermal fits work with global parameters T and fi which, at first sight, appears 
inconsistent with what I just said. Here the role of freeze-out becomes important: freeze- 
out cuts off the hydrodynamical evolution of the thermalized region via a kinetic freeze- 
out criterium |IJ which involves the particle densities, cross sections and the expansion 
rate. In practice freeze-out may, but need not occur at nearly the same temperature 
everywhere M. 

Clearly a thermal fit to hadron production data (if it works) is not the end, 
but rather the beginning of our understanding. One must still check the dynamical 
consistency of the fit parameters Tf, fif, vf. can one find equations of state and initial 
conditions which yield such freeze-out parameters? Which dynamical models can be 
excluded? 

3. The hadronic phase diagram 

In figure [l] I show a recent version of the phase diagram for strongly interacting matter, 
with various sets of data points included 0, || |7], || §]. In the present section I discuss 
the meaning of this figure and explain how these data points were obtained. In the next 
section I will discuss some problems connected with the extraction procedures. 

3.1. Chemical freeze- out 

3.1.1. e + e~ and pp collisions. Let me begin with the e + e~ data point in figure 0. 
(There is also a pp point from Ref. which was omitted for clarity.) In spite of 
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Figure 1. Compilation of freezeout points from SIS to SPS energies. Filled symbols: 
chemical freeze-out points from hadron abundances. Open symbols: thermal freeze- 
out points from momentum spectra and two-particle correlations. For each system, 
chemical and thermal freeze-out were assumed to occur at the same value (1b/T. The 
shaded region indicates the parameter range of the expected transition to a QGP. 
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what I said about case 1. above, a "thermal" analysis of hadron yields in elementary 
collisions [l], |2|] is still interesting. (Of course, in this case the canonical formalism 
must be used, due to the small collision volume.) The interest arises a posteriori from 
the observed universality of the fit parameters, namely a universal "hadronization" or 
"chemical freeze-out" temperature T c h em = ^had ~ 170 MeV (numerically equal to 
the old Hagedorn temperature Th and consistent with the inverse slope parameter 
of the m^-spectra in pp collisions ||), and a universal strangeness fraction A s = 
2(ss)/((uu) + (dd)) (produced ~ 0.2—0.25, almost independent of <Js fl], |2], [10[ . 



This is most easily understood |2| in terms of a universal critical energy density 
e cr it for hadronization which, via the Maximum Entropy Principle, is parametrized 
by a universal "hadronization temperature" Th a( j and which, according to Hagedorn, 
characterizes the upper limit of hadronic phase-space. Supporting evidence comes from 
the observed increase with of the fitted fireball volume Vf (which accomodates 
the increasing multiplicities and widths of the rapidity distributions). Although 
higher collision energies result in larger initial energy densities eo, the collision zone 
subsequently undergoes more (mostly longitudinal and not necessarily hydrodynamical) 
expansion until e cr i t is reached and hadron formation can proceed. The systematics 
of the data can only be understood if hadron formation at e>e cri t (i-e. T>T H for 
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the corresponding Lagrange multipliers) is impossible. With this interpretation, the 
chemical analysis of e + e~, pp and pp collisions does provide one point in the T-^ 
phase diagram (see figure |l|). - The only "childhood memory" of the collision system is 
reflected in the low value of X s , indicating suppressed strange quark production (relative 
to u and d quarks) in the early pre-hadronic stages of the collision. 



3.1.2. AA collisions: strangeness enhancement. In this light the observation || of a 
chemical freeze-out temperature T c h em ~ Th ~ 170 MeV in heavy- ion collisions with 
sulphur beams at the SPS (figure |l|), taken by itself, is not really interesting. It suggests 
that in heavy-ion collisions hadronization occurs via the same statistical hadronic phase- 
space occupation process as in pp collisions. What is interesting, however, is the 
observation [10, 11] that the global strangeness fraction A s ~0.4— 0.45 in AA collisions is 
about a factor 2 larger than in e + e~ and pp collisions. If pp and S+A collisions hadronize 
via the same mechanism, and in S+A collisions the Maximum Entropy particle yields 
fixed at Thad are not modified by inelastic hadronic final state rescattering, this increase 
in X s must reflect a difference in the properties of the prehadronic state! In nuclear 
collisions the prehadronic stage allows for more strangeness production, most likely due 
to a longer lifetime before hadronization. 



It was noted before JTO, O that the global strangeness enhancement occurs 



already in collisions between medium size nuclei (S+S) and remains roughly unchanged 
in Pb+Pb collisions. In this conference we saw data from the WA97 collaboration 
13 , |i~4f which provide two important further details: 



1. While the bulk of the strangeness enhancement from p+Pb to Pb+Pb collisions 
is carried by the kaons and hyperons (A, E), which are enhanced by about a factor 3 
near midrapidity, the enhancement is much stronger for the doubly and triply strange 
baryons 5 and Q and their antiparticles, with an enhancement factor of about 17 (!) 
for Q + at midrapidity. The enhancement clearly scales with the strangeness content 
|i~3|l , as naively expected in statistical and thermal models, but in stark contradiction 
to expectations based on the consideration of the respective production thresholds in 
hadronic (re) interactions. 

2. In semi central Pb+Pb collisions the enhancement grows linearly with the number 
of participating nucleons in the collision, so the enhanced yield of all measured strange 
hadron species per participating nucleon is independent of the effective size of the 
colliding system from about 150 to 400 participants [ P2"f . Where comparison is possible, 
this systematics even carries over to central S+S collisions [IT4J with as few as 55-60 
participating nucleons. So whatever causes the enhancement in Pb+Pb collisions (e.g. 
the existence of a color-deconfined prehadronic stage) is not particular to (semi) cen- 
tral Pb+Pb collisions, but exists already in S+S collisions! I will return to the A- 
independence of this effect below. 

At most half of the 100% increase of global strangeness production in AA collisions 
can be explained |15|] by the removal of canonical constraints in the small e + e~ and 
pp collision volumes (which would be an interesting observation in itself because it 
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would already imply that in nuclear collisions hadron production and the conservation 
of quantum numbers occurs no longer on nucleonic, but indeed on nuclear length scales). 
The remainder of the increase must be due to extra strangeness production in the whole 
fireball volume before hadronization. It is interesting to analyze in the same way the 
strong Q+Cl enhancement in Pb+Pb (by a factor 17 relative to p+Pb JTB)): of course, the 
Q (carrying 3 units of strangeness) suffers a particularly strong canonical suppression 
due to exact strangeness conservation in the small hadronization volume of a pp (or 
p+Pb) collision; for Th a d — 170 MeV, 7 S = 0.5 and V = 17.6 fm 3 (as obtained by fitting 
pp data at y/s = 27 GeV H) Q + Cl are canonically suppressed by a factor 12 relative 
to a grand canonical treatment []lB|. Again the observed enhancement effect in Pb+Pb 
collisions is considerably larger than expected from a simple removal of the canonical 
constraints. 

The observed strangeness fraction X s = 0.45 in nuclear collisions corresponds to 



a strangeness saturation coefficient []T7| 7 S pa 0.7 ||10|| . On the other hand, a value of 
7 S pa 0.7 in the hadronic final state may, in fact, be the upper limit reachable in heavy-ion 
collisions [15| because the corresponding strangeness fraction agrees with that in a fully 
equilibrated QGP at Th a d ~ 170 MeV. If both strangeness and entropy are conserved 
or increase similarly during hadronization, 7 S pa 0.7 in the Maximum Entropy particle 
yield after hadronization would be a universal consequence of a fully thermally and 
chemically equilibrated QGP before hadronization JOJ. The SPS data would then be 
completely consistent with such a prehadronic state. 

The existence of a prehadronic stage without color confinement, both in S+S and 
Pb+Pb collisions at the SPS, is also suggested by an analysis of baryon/antibaryon 
ratios of different strangeness. This was stressed at this conference by A. Bialas [[TjJ] 
who redid this analysis with the new data following the ideas of Rafelski | ]17| |. 

According to figure [l] chemical freeze-out in sulphur-induced collisions at the SPS 
appears to occur right at the critical line, i.e. immediately after hadronization. The SIS 
data, on the other hand, indicate much lower chemical freeze-out temperatures. The 
origin of this is probably due to longer lifetimes of the reaction zone at lower beam 
energies, allowing for some chemical re-equilibration by inelastic hadronic reactions. 

A tendency for some chemical re-equilibration after the hadronization of the 
proposed pre-hadronic stage may also be visible in the still preliminary Pb+Pb data 
at the SPS: although thermal model analyses of these data still give wildly scattering 
results |L0], [7|, [I9[] , some authors |19j find chemical freeze-out temperatures in Pb+Pb 
below 140 MeV. A thermal model analysis of RQMD simulations also gives chemical 
freeze-out temperatures of 172 MeV in S+S, but of only 155 MeV in Pb+Pb collisions 
[2(J. Both analyses show a characteristic failure to reproduce the Q and Cl yields ||12j| . 
This was interpreted |2Tj] in terms of early freeze-out of these triply strange baryons 
due to their small interaction cross sections with other types of hadrons. It is interesting 
to observe that in the thermal analysis of the data JTPJ the model underpredicts the 
measured Q and fl yields (which prefer a higher freeze-out temperature T > 170 MeV 
|0. |l0l) while in RQMD, which is known to produce too few Q and Q baryons in pp and 
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pA collisions, the thermal model overpredicts the simulated yields. All this points to the 
Q and as relatively early hadronic messengers, and the message they seem to carry 
in the Pb+Pb data is again that of the existence of a prehadronic stage with enhanced 
global strangeness which hadronizes statistically at T had = T H — 170 MeV. 



3.2. Flow and thermal freeze- out 

The other important observation in the hadronic sector of nuclear collisions is that of 
collective flow (radial expansion flow, directed and elliptical flow). It is usually extracted 
from the shape of the single-particle momentum distributions. Radial flow, for example, 
leads to a flattening of the mj_-spectra. For the analysis one must distinguish two 
domains. In the relativistic domain p_i_3>mo the inverse slope T app of all particle species 



is the same and given by the blueshift formula ||] T app =TfJ(l+(v±))/(l—(v±)). This 
formula does not allow to disentangle the average radial flow velocity (v±) and freeze- 
out temperature Tf. In the non-relativistic domain piCmo the inverse slope is given 
approximately by T app =Tf+m (v\) , and the rest mass dependence of the "apparent 
temperature" (inverse slope) allows to determine Tf and (vj_) separately, (liapp collisions 
no m -dependence of T app is seen P2"fl .) Plots of T app against m were shown in several 
talks at this conference, showing that the data follow very nicely this systematics, from 
SIS to SPS energies. 

A notable exception are the f2-spectra of WA97 |12] which are steeper than expected 



from this formula. Again, as in the above discussion of their abundance, this reflects 
their character as "early hadronic messengers" [^l|]: the Q and Cl are the only baryons 
which (due to quantum number mismatch) do not have a strong resonance with pions, 
the most abundant particles in the fireball. Since resonance scattering is the most 
efficient thermalization mechanism in a dense hadronic system, the Q and Cl momentum 
distributions freeze out earlier than those of all other baryons. This implies [[H]] that 
they cannot efficiently pick up the collective transverse flow which builds up among the 
pions in the later stages of the expansion, and their spectra reflect the much weaker 
collective transverse flow in the early collision stages, just after hadron formation. 

[This also illustrates the important role of the baryon contamination in the hot 
fireball: in Pb+Pb collisions at the SPS the pion and baryon spectra decouple late and 
cool down to rather low temperatures of about 120 MeV (see below) because the pions 
are "glued" together by the baryons via resonance scattering. At RHIC this glue will be 
less efficient since near midrapidity there will essentially exist only baryon-antibaryon 
pairs at rather low thermal equilibrium abundances. It is thus expected that at RHIC 
thermal decoupling occurs at somewhat higher freeze-out temperatures, closer to the 
hadronization phase transition.] 

The separation of collective flow and random thermal motion from an analysis of 
single particle spectra is not uncontroversial. The main reason is that the fitted values 
for T and v± tend to be strongly correlated. To break the correlation one must study 
spectra of hadrons with different masses in the low-p^ region which, on the other hand, 
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Figure 2. Thermal freeze-out temperature and transverse flow velocities extracted 
from fits to the transverse momentum spectra of negative hadrons (h~) and deuterons 
(d) and to the transverse HBT radius (2tt-BE). The shaded area indicates the overlap 
region near T f . Q . « 120 MeV and v± « 0.55 c. (From Ref. Q.) 




is contaminated by post-freeze-out resonance decays. As a consequence, fits done in 
different p± windows tend to give different results. 

A clearer determination of the transverse collective flow comes from a direct 
measurement of the flow-induced space-momentum correlations via the Mj_-dependence 
of the two-pion HBT radii |'23|| . As shown in [[23], the correlations between 
temperature and transverse flow in a fit of the single particle transverse momentum 
spectra and of the transverse two-particle HBT radii are essentially orthogonal to each 
other (see figure ^), and the combined analysis of spectra and correlations allows for a 
clean separation of random thermal motion from collective flow. 

This analysis [24], ^] gave rise to the open circle for the SPS data in figure 
|l|, indicating the point of thermal decoupling in 158 A GeV/c Pb+Pb collisions 
(Ttherm ~ 120 MeV, (v±) ~ 0.5 c). It is consistent with a comparison of the spectral 



slopes for different mass hadrons by Kampfer [2f|. Earlier analyses of S+S data 
showed a thermal decoupling at T therm « 140 — 150 MeV, (v±) ~ 0.25 — 0.35 c. The open 
circle for the AGS as well as the open quadrangles for the SIS in figure [I] were obtained 
similarly as in [[26], 2^] by comparing transverse momentum spectra of particles with 
different masses (see |7|, [8]] for references). All open symbols correspond to heavy collision 
systems (Pb+Pb, Au+Au). The dashed line connects them by eye in an attempt to 
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construct a "thermal freeze-out curve" for heavy-ion collisions of size 200+200. 

It is interesting to analyze the A-dependence of the various hadronic observables, 
i.e. of strangeness enhancement, chemical and thermal freeze-out temperatures and 
radial transverse flow. This discussion should include the available information on the 
impact parameter dependence in Pb+Pb collisions since collisions at different impact 
parameters also involve different numbers N paxt of participating nucleons. Whereas the 
freeze-out temperatures (certainly the thermal freeze-out temperature, but perhaps also 
the chemical one, see above) seem to come down with increasing size of the collision 
system, while at the same time the strength of the radial transverse flow goes up, the 
strangeness enhancement (i.e the strange particle production per participating nucleon 
relative to pp and pA collisions) appears to be independent of the number of iV part , at 
least in the range 60 < iV par t < 400. The buildup of radial flow is largely (although not 
exclusively) a hadronic reinteraction phenomenon |]2T|; the same is true for the freeze- 
out process. The iV part -dependence of both features can be explained in terms of the 
longer lifetime of the reaction zone as N part increases, giving the system more time to 
re-equilibrate, cool down and develop collectivity. 

In contrast, the A-independence of the strangeness enhancement features suggest 
that they are not due to hadronic re-interactions, but originate in a prehadronic phase 
with properties which are essentially independent of the system size once iV pai . t > 50 or 
so. My interpretation of these facts is that at SPS energies the energy density threshold 
for QGP has been overcome by a sizeable margin f2l| , and that even in small collision 
systems a deconfmed phase is created which interacts sufficiently long and sufficiently 
strongly to approximately saturate strangeness production. One can even argue that 
isotropic pressure (a signature of local thermalization) must be present at this early stage 
P8fl , and that the observed elliptic flow in non-central Pb+Pb collisions |2P| actually 
signals this early pressure f30fl . A future systematic investigation also of the range 
1 < iVpart < 100 would be very useful to study the onset and saturation of thermal and 
collective behaviour as the size and lifetime of the collision system increases. It will not, 
however, be an efficient method to study the onset of deconfinement - for that one must 
go to lower beam energies. 



4. Limitations of thermal model analyses 

After having explained how the data points in figure |l| were obtained, I would now like to 
ask the notorious David Mermin question "What's wrong with this phase diagram?" 
In other words, I want to point out in more detail certain unavoidable problems 
with thermal model analyses of heavy-ion data. Only by remaining conscious of the 
limitations of the thermal approach and avoiding the overinterpretation of uncontrollable 
details one can fully exploit its power in providing essential qualitative insight into the 
physics of heavy-ion collisions. 
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4-1. Rapidity dependence of particle ratios 

The first problem is of purely practical nature: no heavy-ion experiment so far has full 
Arc acceptance for identified particles, and particle spectra are available in restricted 
windows of p± and y only. The observed nearly exponential form of the p_i_-spectra 
allows for an extrapolation of the yields to the full transverse momentum range without 
introducing large uncertainties (at least if the acceptance covers sufficiently low values 
of p±). A similar extrapolation of the rapidity spectra is not possible: even in a system 
which is in perfect local thermal equilibrium, particles with different masses tend to have 
strongly different rapidity distributions. In practice these differences in the shape of the 
rapidity spectra are even stronger, in the sense that even particles and antiparticles 
(with obviously identical masses) have different rapidity distributions. Thus there is 
essentially no way of extrapolating the rapidity distributions without really measuring 
them. 

That this is a serious problem for thermal model analyses is illustrated by the 
following example: consider a stationary, spherical fireball in global thermodynamic 
equilibrium. It emits hadrons with the following rapidity distributions: 

x 2" 



^ e - mi co8hy/T 



dy 



T T 
1 + 2 — + 2 



m i cosh y V m j cosh y 



(1) 



These resemble Gaussians with widths Tj m 2.35y^- (the approximation being valid 
for nii 3> T). Clearly, the particle ratios (dNi/dy) / (dNj /dy) then depend strongly on 
the position of the rapidity interval dy: away from y = heavy particles will be much 
more strongly suppressed relative to light particles than near y — 0. In this case a 
measurement of particle yields in a small rapidity window is completely useless for a 
thermal model analysis (irrespective of whether the window is located near y = or at 
y 7^ 0) unless it is known a priori that the radiator is a stationary spherical fireball. 

The presence of (strong) longitudinal flow in relativistic heavy-ion collisions does 
not help very much in this connection; only in the limit of infinite beam energy with 
exact longitudinal boost-invariance due to Bjorken scaling, resulting in flat rapidity 
distributions, is it possible to base a thermal analysis on data in a finite, narrow rapidity 
interval. At SPS energies and below, where Bjorken scaling is not observed, thermal 
fits of data in finite rapidity windows require that the thermal model yields are cut to 
the actual experimental acceptance; this induces serious dependences on the detailed 
model assumptions, for example about the strength and profile of the longitudinal and 
transverse flow of the source. This mo del- dependence was recently studied in some 



detail in Ref. 32 



Flow effects drop largely out, however, if one works with particle ratios obtained 
from Air yields. (Please note that this requires a measurement of some sort, not a blind 
extrapolation of data from a small window in y and p± to full momentum space!) The 
insensitivity to hydrodynamic flow becomes exact if the freeze-out temperature and 
chemical potential is everywhere the same. If freeze-out occurs on a sharp hypersurface 
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E, the total yield of particle species i is then given by 

Ni = Ijf**^ f ^ p) = jj 3 ^)jt(x) , (2) 

where 

jjftx) = J <fip20(pW - m*)p ^ Mx) 9 ^ /T±i (3) 

is the number current density of particle species i. In thermal equilibrium it is given by 
ji(x) = p i (x)u fl (x) with 

Pi{x) = u^(x)j^(x) = J d 4 p26(p°)5(p 2 - m*)p-u(x) fi(p-u(x);T, p^ 

= J d 3 p' f t {E p ,- T, /i,) = Pi(T, p^. (4) 

Here E p > is the energy in the local rest frame at point x. The total particle yield of 
species i is therefore 

N i = p i (T,fH)J <Pa^x)u*(x) = p i (T,to)Vv(u' t ) (5) 

where only the total comoving volume of the freeze-out hypersurface E depends on 
the flow profile u^. Thus the flow pattern drops out from ratios of Air yields which 
therefore depend only on T and the chemical potentials. These considerations are easily 
generalized to "fuzzy freeze-out" (i.e. freeze-out from a space-time volume rather than 
from a sharp hypersurface): as long as T and pi are the same everywhere, 4ir particle 
ratios are independent of the collective dynamics of the source. 

For heavy-ion collisions at SPS energies and below one should therefore perform a 
thermal analysis on 47r-integrated yields and not on particle ratios inside small rapidity 
windows. This requires a strong experimental effort to measure the rapidity distributions 
of as many particle species as possible over the full rapidity range. 



4-2. N on- constant thermodynamic parameters at freeze- out 



The second, even more serious problem is the observation that in reality freeze-out does 
not happen at constant temperature and chemical potential. For example, it was shown 
in Ref. |33j that a successful thermal description of the rapidity distributions of hadrons 
created in 200 A GeV S+S collisions (in particular the different shapes of the rapidity 
distributions for A and A, K + and K~) requires not only strong longitudinal flow, but 
also a baryon chemical potential Pi{rf) which depends on the longitudinal position in 
the fireball: the central rapidity region is baryon-poorer than the target and projectile 
fragmentation regions. A second example demonstrating this type of problem is the 
observed rapidity dependence of the p^-slopes of the h~ and proton spectra |34j and of 
the ii^-slopes of the transverse HBT radius R± [|4], |3~5|| . According to a simultaneous 
analysis (as discussed in section |37^ above) of spectra and correlations from Pb+Pb 
collisions at the SPS by Schonfelder [35| the decrease of the inverse slope parameters 



away from midrapidity must be attributed to both a reduction of the transverse collective 
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flow and of the freeze-out temperature T{rj) . If true, this would speak against a constant 
freeze-out temperature. 

In such a situation the thermal fit replaces the functions T(rj) and by suitably 
averaged values T, fii (see Ref. |32| for a recent detailed study). Obviously the fit 



will then not be perfect: particle yields from a system in perfect local thermodynamic 
equilibrium (as, e.g., assumed in all hydrodynamic simulations), but with spatially 
varying temperature and chemical potentials, cannot be exactly recovered by a fit with 
constant temperature and chemical potential. In practice, this does not appear to be a 
very serious problem if one believes that the results from a recent thermal model analysis 
of particle yields from a hydrodynamic simulation [3^] are representative for realistic 



situations: the freeze-out temperature reconstructed from the fit nearly coincided with 
the input temperature at which the hydrodynamic evolution was stopped, and the fitted 
chemical potentials agreed approximately with their average values across the freeze- 
out surface. Nevertheless, small differences remain between the real yields from the 
hydrodynamic simulation and the yields returned by the thermal model fit. In a least 
mean square fit these systematic deviations would lead to a value for x 2 /d.o.f. which 
increases above all limits as the statistical error of the simulated ("measured") yields 
is further and further reduced, even though the system was, by construction, in perfect 
local thermal equilibrium. 

This illustrates that global thermal fits to heavy-ion data can never be fully 
successful, due to the dynamics of the collision and its intricate influence on the freeze- 
out process. For this reason one must not expect too much from the thermal model - 
a perfect fit with extremely small x 2 /d.o.f. is not necessarily a good and often rather 
a bad sign, indicating accidental error correlations, e.g. due to the use of redundant fit 
parameters. 

This leaves us with the question where to draw the line between "good" and "bad" 
thermal model fits, between success and failure of a thermodynamic description of 
relativistic heavy-ion collisions. The above discussion should have made it clear that 
X 2 /d.o.f. is not a good criterium for answering this question. On the other hand, 
a fit which reproduces particle yields which cover a range of more than three orders 
of magnitude with individual deviations of less than ±25% @ is obviously not bad. 



Quantitative model studies like those presented in |32| give us guidance for separating 
the grain from the straw; when supplemented by a thermal analysis of microscopic 
kinetic simulations as those presented at this meeting by Larissa Bravina |36| , they are 
the foundations which we can use when collecting arguments in favor or against the 
formation of thermalized hot hadronic matter and quark-gluon plasma. 

5. Conclusions 

Let me summarize shortly: a thermal + flow analysis of yields, spectra and 2-particle 
correlations in S+A and Pb+Pb collisions at the CERN SPS suggests 

• the formation of a prehadronic state in which twice as much strangeness is produced 
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as in pp and pA collisions and in which quarks are uncorrected, i.e. they are not 
bound into color singlets; 

• statistical hadronization of this state at Th a( j ~ 170 ± 20 MeV with hadron 
abundances controlled by the Maximum Entropy Principle; 

• rapid decoupling of the particle abundances, with chemical freeze-out temperatures 
Tchem ~ Thad m sulphur-induced and T c h em < ^had in Pb+Pb collisions; 

• elastic rescattering among the hadrons (dominated by s-channel resonances) after 
hadronization from the prehadronic state which leads to a cooling of the momentum 
spectra and the generation of (more) collective flow; 

• finally, thermal freeze-out at T t herm ~ 140 — 150 MeV in S+A and at T t herm ~ 
120 ±10 MeV in Pb+Pb collisions, with average transverse collective flow velocities 
of order (v±) ~ 0.4 — 0.5 c. 

Smaller collision systems distinguish themselves from larger ones not primarily by 
the achieved maximal energy density, but by the occupied collision volume in space and 
time. Compared to S+S collisions, Pb+Pb collisions live longer until thermal freeze-out, 
expand more in the tranverse direction, develop more transverse collective flow and cool 
down to lower (chemical and thermal) freeze-out temperatures. Of course, this view 
disagrees with the (present) majority opinion (I refer to the respective contributions 
to the Proceedings of "Quark Matter '97" that the critical energy density for 

deconfinement can be crossed at fixed beam energy of 160 A GeV by changing the size 
of the projectile and target or the collision centrality, and that the "anomalous" J/i/j 
suppression observed in central Pb+Pb collisions as a function of produced transverse 
energy |37], [$8| signals this transition. 



I find this interpretation irreconcilable with the systematics of light hadron 
production as discussed in this talk; a more consistent interpretation rests on the 
observation |38| that, as one increases iVp art , one first sees "anomalous" suppression 
of the weakly bound xp', then (in semiperipheral Pb+Pb collisions) the "anomalous" 
suppression of the more strongly bound Xc states (indirectly, via the disappearance of 
their 32% feed-down contribution to the measured J/ip yield), and only in very central 
Pb+Pb collisions the disappearance of directly produced J/ip's which are very strongly 
bound (this last part of the supression pattern still remains to be confirmed by an 
improved measurement at very high Et)- This suggests to me that what NA50 is 
seeing is not the onset of deconfinement (the latter is there even in S+A collisions), but 
the dissociation of more and more strongly bound heavy quark states (respectively the 
removal of the corresponding components in the cc wavefunction) by collisions with the 
dense partonic medium in the early stages of the collision. For the more strongly bound 
states most of these collisions will be subthreshold; for this reason a longer lifetime of the 
dense early stage, which is achieved in larger collision systems or more central collisions, 
is crucial for an efficient destruction not only of the weakly bound i/j', but also of the 
more strongly bound Xc and J/ip. Charmonium suppression is thus, in my opinion, 
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more a lifetime effect than a deconfinement signal (although the necessary high density 
of scatterers with sufficiently large cross sections probably requires deconfinement, too). 

The picture which I have painted is, I believe, intrinsically consistent. It is 
sufficiently simple to be attractive but also sufficiently sophisticated not to be unrealistic. 
It may not be unique, not least because of the intrinsic systematic uncertainties 
associated with thermal model analyses which I pointed out and not all of which are 
quantitatively understood. What is urgently needed is more high-quality data on the 
chemistry of Pb+Pb collisions, the reconciliation of some puzzling discrepancies between 
different experiments as discussed at this meeting, and an improved systematics of the 
impact parameter and A-dependences, both in the light and heavy hadron sector. 
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